An unification of general theory of relativity with Dirac's large 

number hypothesis 

H. W. Peng 1 ' 

1 Institute of Theoretical Physics, Chinese Academy of Sciences, 
P.O.Box 2735 Beijing 100080, China 

Abstract 

Taking a hint from Dirac's large number hypothesis, we note the existence of cosmologically 
combined conservation laws that work to cosmologically long time. We thus modify Einstein's 
theory of general relativity with fixed gravitation constant G to a theory for varying G, with a 
tensor term arising naturally from the derivatives of G in place of the cosmological constant term 
usually introduced ad hoc. The modified theory, when applied to cosmology, is consistent with 
Dirac's large number hypothesis, and gives a theoretical Hubble's relation not contradicting the 
observational data. For phenomena of duration and distance short compared with that of the 
universe, our theory reduces to Einstein's theory with G being constant outside the gravitating 
matter, and thus also passes the crucial tests of Einstein's theory. 
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In Einstein's theory of general relativity as in Newton's theory of gravitation, the 
strength of the gravitational interaction is described by a fixed dimensional constant Gn = 
6.7 x 10- 11 m 3 kg- 1 s -2 = 7.4 x KT^mkg- 1 = 8.2 x KT^mJ- 1 by taking c = 3 x 10 8 m/s= 1. 
Einstein's theory has been applied to cosmology; another cosmological constant with di- 
mension of length to the power minus two was once specially introduced by Einstein in 
his attempt to construct a static model for the universe. This soon has lost his favor as 
observations of the red shifts of extra-galactic nebulae show definitely that the universe is 
expanding, but cosmological constant is still used today as a free parameter in trying to 
fit the observational data with the non-static homogeneous cosmological model using the 
Roberston- Walker metric. Rough estimates have thus been obtained about seventy years ago 
for the age t of universe at present to be of the order of 10 9 years, and recently of the order 
of 10 10 years. This is about 10 39 or 10 40 times the time (e 2 /m e c 3 ) needed for light to travel a 
distance of the classical radius of the electron. Similarly, with Hubble's estimate p = (1.3 to 
1.6) x 1CT 30 g/cm 3 for the density of matter due to the extra-galactic nebulae averaged over 
cosmic space and a factor thousand or hundred times to include various invisible matter, one 
can estimate the ratio M/m p of the mass M of matter in the universe of radius R at present 
to the proton mass m p , and obtain another large number of the order 10 78 or 10 80 which 
is about the square of the pervious large number. Further, the ratio e 2 /(Gm p m e ) of the 
electrostatic to the gravitational force between the proton and the electron in an hydrogen 
atom is also a large dimensionless number of the order of 10 39 . On comparing these num- 
bers, Dirac proposes Q the large number hypothesis - a sort of a general principle that very 
large numbers which turn up in Nature and have no dimensions are related to each other. 
Accepting this hypothesis and comparing the above three large dimensionless numbers we 
obtain the following two large number equations where small number coefficients are omitted 



Dirac concludes [l| from ((TJ, (J2J) that G oc t^ 1 and Mat 2 during the evolution of the 
universe, assuming m e , m p and e be constant. Then neither M nor G separately but the 
combination G 2 M conserves in cosmologically long time. If we take (|T|). as empirically 
true, we have on dividing (j2j) by ((TJ a relation free from e, m e , m p , 




(1) 
(2) 



M/m p 



{t/(e 2 /m e c 3 )} 2 . 



GM « cH « t, 
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on putting c = 1. This exhibits clearly the inconsistency of applying Einstein's theory to 
cosmology, because in Einstein's theory G is constant and mass M conserves, but the age 
of universe t does vary. 

In this Letter, we take the existence of the cosmologically combined conservation law as 
a hint in our attempts to generalize Einstein's theory to a theory with varying G, working 
consistently with the above large number relations to cosmologically long time. In contrast 
to Dirac's assumption of constant m e , m p and e , we believe in the unity of physical laws, 
so that m e , m p will evolve in the same way as M, i.e. Gm oc t by Q and hence e 2 oc tm e by 
(J2J). Instead of Dirac's conclusion G oc t -1 , M oc t 2 we can only conclude that if we assume 
that G oc t~ n , then m oc t 1+n , (m = m e , m p or M), and e oc f 1 + n / 2 with n likely but not 
necessarily equal to one. If we make use of the constants t and G = G(t ), we can write 
the proportion as an equality. Introducing a new dimensionless variable c/> 2 defined by the 
n-th root of G/Gq so that our fundamental assumption becomes 



The above proportion laws of evolution are then converted into cosmologically combined 
(cc.) conservation laws that do remain constant over cosmologically long time. 



We shall call such quantity defined by multiplication with an appropriate power of </> the 
cosmologically combined (cc.) quantity; it is of the same dimension as the original quantity, 
and its numerical value is equal to that of the latter at the epoch to, since the value of c/> 
at the epoch t is identically equal to one. If we choose t to be the epoch at present, then 
we may use Gn for G$. We shall denote the cc quantity with a decoration on the letter as 
exemplified in (j^J). 

We modify Einstein's theory of general relativity (here we follow the treatment by Dirac 
P] with comprehensive action principle) as follows. For matter we change the action integral 
into 




(4) 




(5) 




(6) 



with the modified constraint in the form of an ordinary divergence relation 




(7) 
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where (we use the letter u for the four- velocity instead of Dirac's letter v) 



j? = u^(-g), (8) 

with 

P = 2 " + V (9) 
This constraint gives the conservation of the c.c. mass m for homogeneous <fi = 4>(t) 

m = I p i dx 1 dx 2 dx 3 = <p 2n+2 m. (10) 



Thus our modification of Einstein's theory consists of inserting a factor </> 2n+2 in the action 
integral and introducing c.c. variables. Since Einstein's original gravitational action integral 
contains a denominator G = Go<f) 2n , so we take for the modified gravitational action integral 
the expression where only a factor 2 is left after the cancellation 

I g = (16-kGo)- 1 J [4?R pv - wfrMg^i-g^x. (11) 

Here we have added a kinetic term to account for the <p to vary, the form of the kinetic term 
being unique because forms like 4>4> w , v can be transformed into ;At ;t / by an integration by 
parts. The numerical coefficient w of this kinetic term will be chosen to be w = 8, as will 
be shown below to correspond with the case of k = for the Robert son- Walker metric for 
homogeneous cosmological models. We note that our expression pijl . though it looks like 
that which has been used in Brans-Dicke theory of varying G, is actually different from the 
latter because of the different relation between and G by our having inserted the factor 
(f) 2n+2 in the action integral. The comprehensive action principle 5I to t = SI g + 5I m = gives 
the variational equations for g^ u and (f), 

= ^{K - (i/2)i#$a + mi - (^y-x - m<m ;/3 - (i/2)</v<r^] 

= -87tG tJ = -%nG Q pu a vP, (12) 
$ = (W£ + ^0) = O, (13) 



and the equations of motion pu v u^ v = together with the constraint (p-u M ) ;M = for the 

i _ 



elements of matter. These equations are not all independent. Owing to the identities T^. 



0, which follow from the equations of motion and the constraint for matter, and N^.p+t 
0, which can be derived from the invariance of the action integral I g by an infinitesimal 



coordinate transformation as indicated in Dirac's book |2j, or directly verified from ()12|) ()13|) 
by tensor calculus, these equations are compatible but indeterminate to allow to change 
with coordinate transformation. We note that in our theory the following combination of 
(JI2J) and (H3J) 

+ 0$ = (w/2 - 3)(0 2 )£ = -8vrG p (14) 

is particularly simple and can be used as one of the independent equations. Also it is 
interesting to write the field equations for g^ in a form analogous to that in Einstein's 
theory by dividing both sides of (fT2j) by <ft 2 : 

Ri - (l/2)R a X + Ki = -87iG pu a uP/(P 2 = -8nGpu a u^. (15) 

Here on the right G is not a constant but varies as 4> 2n while on the left appears the tensor 
which is completely determined by the derivatives of 4> 

K = [(<I>X - (0 2 );^]/0 2 

-w[<P, a ^-(l/2)<p, a r5§/<p 2 (16) 

and differs in its tensor character from the cosmological term A5@ usually introduced ad 
hoc in Einstein's theory at the same place. Thus in our theory there occurs a variable 
cosmological tensor but no ad hoc cosmological constant A as a free parameter. 

We now apply our theory to homogeneous cosmological model. We can adopt as usual the 
Robert son- Walker metric because in simplifying the metric to this form only considerations 
on symmetry and freedom of coordinate transformation have been used, but no use is made 
of the field equations. We have, as usual 

ds 2 = dt 2 - R 2 (t){dr 2 /(1 - kr 2 ) + r 2 d6 2 + r 2 sin 2 6d<p 2 } (17) 

with t denoting the cosmic time and r, 9, if the dimensionless co-moving coordinates. For 
the Robertson- Walker metric the three non-trivial field equations [one from (|13|) and two 
from (fl~2|) according to a = (3 = 4, and a = [3 ^ 4] are connected by one identity (a = 4), so 
there are only two independent equations to determine the two unknown function <j) 2 (t) and 
R(t). We shall choose (fT2|) with a — (3 — 4 and (fT4]l as the two independent field equations. 
All the equations of motions for the element of matter are trivially satisfied by their being 
at rest in the co-moving coordinates, i.e. u l = Ui = for i = 1, 2, 3 and hence u A = = 1. 
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From the modified constraint we obtain p m = p m (t')R 3 (t')/R 3 where the constant argument 
t' is arbitrary. Then (j!4|) can be integrated to give 

R 3 d((j) 2 )/dt = (3 - w/2)- 1 87vG p m (t')R 3 (t')t, (18) 

the spatial factor involved in y/ (—g) being cancelled out. The constant of integration additive 
to t is chosen to be zero so that t is the age counted since the big bang when R = O.With 
the help of c/> 2 oc £ -1 , we obtain from (fTH)) that R oc t, i.e. the rate dR/dt = R of 
expansion is a constant. The other independent equation, namely (|T2|) with a = (3 = 4, is 
for Robertson- Walker metric explicitly 

2 [-3(i? 2 + k)/R 2 } - 3(<f) 2 )R/R - w<j) 2 /2 
= -8nG p(t')R 3 (t')/R 3 . (19) 

Writing proportions as equalities, for an arbitrary t' 

2 = t'<p 2 (t')/t , R/R(t') = t/t', (20) 

we see that (QSJ is satisfied with 8irG p(t')(t') 2 /(p 2 {t') = (w/2- 3) and (HU is satisfied 
with 8TiG p{t'){t') 2 /<p 2 {t') = w/8 + 3k/(dR/dt) 2 . These conditions agree only when w/8 = 
1 + k/ '(dR/dt) 2 . These conditions also show with and that p(t) oc t~ 2+n so we have 
< n < 2 for p(t) to decrease with t increasing. If we believe that (dR/dt) 2 < 1, the choice 
of w/8 = 1 determining thereby k = seems more reasonable than any other value chosen 
among the range w/8 > 2 for k = +1, while the case of k = — 1 is rejected by (JTSJ). 

With w = 8 and so k = 0, the solution (}2T)j) gives by (fT9"j) = 2/t 2 , and the equation 
(PT)|) can be expressed as Q m + = 1 with Q m = 1/3 and £7a = 2/3 in the notation 
familiar with cosmologists. The other non-varnishing elements of our cosmological tensor 
are Aj = A 2 , = A 3 , = l/t 2 , being one half of A4. 

With k = and R oc t, we can derive the theoretical Hubble relation very simply by 
following Weinberg's (jj treatment. The luminosity distance di is given by the general 
formula di = (1 + z)riR(t ) with z = R(t ) / R(t±) — 1, the suffix referring to the observer 
at the origin tq = 0, the suffix 1 refereing to the source at r\. In our case with r±R(to) = 
R(to) Jq 1 dr = R(t ) f*° dt/R = t log(t /ti) = t log[R(t ) / R(h)], we obtain the theoretical 
Hubble relation 

d L = t (l + z)\og(l + z), (21) 
6 



with the Hubble constant at present Hq = t^ 1 of course. This relation is compared with 
the recent experimental data of Riess et al [4( up to z = 1.55 by Dr. Wang, as shown in 
Fig^ The value of the Hubble constant 61.2 kms~ 1 Mpc~ 1 obtained by means of least x 2 
falls within the range Ho = 64 ± 3 kms~ 1 Mpc~ 1 of their early analysis for z < 0.1 5]. The 
corresponding to = Hq 1 is 16 Gyr which is consistent with the 1 — a fit by Tegmark with 
r < 0.3 fl. 

Turning to the crucial tests of Einstein's theory, which all concern with phenomena 
of short duration and small distance outside the gravitating body, e.g. our Sun, at the 
present epoch t - We shall consider more generally similar phenomena at the epoch t e , 
say. Neglecting the difference t — t e against t e we can use the static or, to be precise, 
the quasi-static solution of (fT2"j) . (fT^j) by neglecting all time derivatives. For our purpose 
we need only the exterior solution with p = on the right hand side of these equations. 
Then from (|14j) that we see that <fi 2 is independent of the spatial variables, with a value 
determined by the boundary condition at far from the gravitating body at the same epoch, 
(j) 2 = <f) 2 (t e ). Then ()12j) or (fTB^l reduces to that in Einstein's theory, as all spatial derivatives 
in the cosmological tensor vanishes for constant and time derivatives are neglected. Thus 
in our theory the static metric outside the gravitating body will be the same as in Einstein's 
theory. If we treat the sun as a mass point at rest with p = mS(x)S(y)5(z), we obtain the 
static exterior Schwarzschild solution where the constant of integration, usually denoted by 
Gm, actually comes from the coefficient of 5 functions in Gom5(x)5(y)5(z)/(j) 2 , so it is the 
quasistatic evolving value G(t e )m(t e ) at the epoch t e , as shown by (HJ), ©• To compare 
with the crucial tests of Einstein's theory, we take the present epoch to for t e and use for 
Gm the correspond value at the present epoch, which is the same value used in Einstein's 
theory. Thus we obtain the same exterior metric in our theory, as in Einstein's theory, 
and quasi-static spatially constant 4>(to) = 1 outside the gravitating body. Hence the three 
crucial tests, the advance of Mercury's perihelion, the deflection of light by the Sun, and 
the gravitational red-shift of the Sun, are all passed by our theory as by Einstein's theory, 
because of the same geodesies, the same null geodesies, and the same g 44 . 

Finally I want to express my indebtedness to Prof. Yu Yunqiang for drawing my attention 
to the luminosity distance, to Dr. Wang Xiulian for making comparison with the data and 
contributing the figure, and to Tu Zhanchun for making this manuscript to revtex4 style. 
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FIG. 1: The luminosity distance and redshift. The experimental data is taken from Ref.jj] and 
the solid line is fitted by Ea.lf2T|). 
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